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INTRODUCTION 


The  molecules  H^O,  NO^,  O^,  CO2  and  N^O  are  the  dominant  triatomic  species 
in  the  earth's  atmosphere.  Because  of  the  ready  availability  of  high  energy 
electromagnetic  radiation  in  the  upper  atmosphere  their  behavior  can  be  quite 
diverse.  Absorption  processes  leave  the  molecules  in  highly  excited  and  ioni- 
zed states,  and  in  these  cases  geometries  often  differ  drastically  from  the 
ground  state  geometries.  Furthermore,  dissociation  and  predissociation  to  atom- 
diatom  asymptotes  is  common.  In  order  to  understand  such  atmospheric  processes, 
certain  fundamental  quantities  are  needed;  among  these  are  the  definitive  loca- 
tion of  the  excited  energy  levels  and  the  delineation  of  reaction  mechanisms. 

These  same  quantities  clearly  are  necessary  for  describing  the  basic  che^'ical 
behavior  of  these  molecules.  In  spite  of  this  need,  most  existing  information 
is  inadequate  for  interpretive  or  predictive  purposes,  and  understanding  atmos- 
pheric reactions  remains  largely  guesswork.  Our  investigation  is  concerned  with 
providing  these  energy  levels  and  reaction  surfaces  through  the  use  of  well  esta- 
blished ab  initio  techniques.  This  new  theoretical  knowledge  will  be  used  to 
advance  experimental  analysis  by  participating  in  a feedback  process  which  histori- 
cally has  proved  to  provide  a deeper  understanding  of  the  phenomena  involved.  We 
hope  that  the  foundation  established  by  this  research  will  lead  to  an  expanded 
effort  in  atmospheric  computations  at  Argonne.  Our  concerns  the  past  year  have 
largely  been  excited  state  characterization  and  method  calibration.  Both  are 
prerequisites  for  any  long-range  investigation.  This  is  particularly  true  of 
the  air  triatomics  because  there  is  so  little  definitive  information  presently 
available.  We  have  also  determined  linear  geometry  to  bent  geometry  correlation 
diagrams.  These  provide  simple  ways  of  predicting  reaction  possibilities.  For 
example,  we  found  that  the  four  lowest  excited  states  of  CO^  are  most  stable  in 
linear  geometries.  The  fifth  excited  state  was  found  to  be  more  stable  in  a bent 
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geometry.  Moreover,  in  the  bent  geometry  it  becomes  more  stable  than  three  of 
the  four  lowest  linear  ion  states,  and  iience,  very  likely  predissociates  the 
lower  energy  linear  states  of  CO^.  It  is  probable  that  the  predissociation 
just  described  for  CO^  carries  over  into  the  Rydberg  series  that  ultimately 
converge  to  the  ion  states  in  question.  Both  types  of  predissociation  could 
be  important  in  atmospheric  processes.  There  are  numerous  similar  examples  of 
important  information  being  provided  by  this  theoretical  study. 

We  have  used  a variety  of  quantum  chemistry  computational  techniques  in 
our  study  of  the  atmospheric  triatomic  molecules.  Certainly,  an  integral  part 
of  our  systematic  study  is  the  calibration  of  various  levels  and  methods  of  com- 
putation. The  computational  methods  used  here  include:  (1)  expansion  or  Ro- 

othan  self-consistent  field  (SCF)^,  (2)  traditional  configuration  interaction 

2 

(Cl)  using  canonical  molecular  orbitals  , (3)  multi-configuration  self-consis- 

3 4 5 

tent  field/configuration  interaction  (MCSCF/CI)  ’ ’ , and  (4)  equations-of- 

fi  7 

motion  (EOM)  . The  computer  cedes  used  in  our  study  were  (1)  the  BISONMC 

program  developed  by  Das  and  Wahl  for  computing  SCF  and  MCSCF  wavefunctions,  (2) 

the  Battel le-Ohio  State  Cl  package  developed  by  Shavitt  and  co-workers,  and  (3) 

the  EOM  programs  developed  by  McKoy  and  co-workers  at  Caltech.  The  following 

subsections  consist  of  a molecule-by-molecule  summary  of  the  results  which  are 

the  body  of  this  report. 

is  the  simplest  prototype  for  ion-molecule  reactions  among  the  air 

2 

triatomics.  We  have  computed  potential  energy  surfaces  for  the  ground  B-j  and 

2 + 2 
the  B2  states  of  H2O  . The  results  of  our  vibrational  analysis  of  the  B.j  sur- 

2 

face  are  in  good  agreement  with  experiment.  The  Bo  surface  is  thought  to  be 
predissociated  twice,  and  has  a complicated  anharmonic  vibrational  structure. 

The  present  level  of  our  results  support  either  of  two  experimental  claims  re- 
garding the  band  origins. 
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COq.  Previous  theoretical  spectra  were  shown  to  be  of  limited  accuracy.  Experi- 
mental assignments  based  on  these  are  therefore  not  definitive.  Our  best  corre- 
lated spectra  are  the  most  accurate  available  for  the  lowest  multiplet  of  CO^. 

Linear  geometry  to  bent  geometry  correlation  diagrams  were  also  determined  for 
the  CO2  and  CO2  systems.  Excited  states  of  CO^  were  found  to  be  stable  in 
bent  geometries  relative  to  atom-diatom  dissociation.  Several  of  the  excited 
linear  states  in  CO^  are  likely  to  predissociate  via  a lower-lying  bent  geometry 
state.  Conflicts  with  existing  CO2  level  assignments  were  found,  and  the  spectrum 
was  reassigned.  Finally,  the  character  of  vertically  excited  states  in  neutral 
CO2  was  determined. 

O^.  The  ongoing  controversy  over  the  reaction  of  ozone  with  photolyzed  freons 
underlines  the  important  role  these  triatomics  play  in  atmospheric  chemistry. 

Our  correlated  results  for  the  vertical  excitation  energies  of  0^  are  in  very 
good  agreement  with  experiment,  and  with  other  calculations.  The  0^  and  0^  cal- 
culations provide  qualitatively  correct  electronic  spectra  for  the  ions.  The 

ordering  of  states  in  0^  is  controversial  at  present,  and  we  support  the  ordering 
2 2 2 

X Ap  A 82  and  B A2.  Our  0^  spectrum  is  the  only  existing  data  concerning  the 
excitation  energies  of  this  ion. 

f^20-  The  ground  state  potential  energy  surface  of  N20~  is  characterized  and 

related  to  the  ground  state  potential  energy  surface  of  the  neutral  N2O  molecule 
by  combining  theoretical  and  experimental  information.  Among  the  quantities  re- 
ported are  equilibrium  geometries,  dissociation  energies,  vertical  and  adiabatic 
electron  affinities,  and  the  minimum  intersection  locus  of  the  ground  state  N2O  j 

and  N2O”  surfaces.  The  results  imply  that  the  reaction  0”  + N2  *■  N2O  + e will 
be  strongly  enhanced  by  N2  vibrational  excitation. 

NO 2.  A definitive  theoretical  study  of  N02^was  performed  earlier  in  our  labora- 
tory. Agreement  with  known  experimental  results  was  very  good.  The  energy  levels 
obtained  in  our  present  effort  provide  similar  data  for  NO^  and  NO^,  where  little 

_ J 
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reliable  information  is  available.  Agreement  with  our  previous  NO^  study  is 
good.  We  suggest  that  known  spectral  complexities  in  NO^  are  due  to  vibration- 
ally  excited  states.  We  also  found  that  existing  theoretical  NO^  spectra  do  not 
consider  all  of  the  electronic  states  in  their  reported  energy  ranges,  and  hence, 
spectral  assignments  based  on  these  are  unreliable.  Our  NO^  energy  levels  are 
the  only  existing  ah  initio  information. 

Part  of  our  computational  efforts  were  supported  by  the  CDC6600  computers 
at  Wright-Patterson  Air  Force  Fase.  We  discuss  our  experience  performing  large 
scale  molecular  calculations  remotely  on  the  Wright-Patterson  computers  in  Ref- 
erence 8 of  the  publication  list  given  in  the  next  section.  Many  elements  of 
the  effort  are  general  and  can  be  expected  to  be  encountered  by  other  researchers. 
A summary  of  the  historically  important  role  played  by  Wright-Patterson  is  also 
given. 

Personnel  Changes.  At  the  end  of  the  last  fiscal  year.  Dr.  Richard  C.  Raffenetti 
replaced  Dr.  P.  J.  Fortune  (who  joined  the  staff  of  the  Applied  Mathematics  Divi- 
sion at  Argonne)  as  project  scientist  on  the  Atmospheric  Triatomic  Molecules 
project.  Dr.  Raffenetti  received  his  Ph.D.  in  physical  chemistry  from  Iowa 
State  University  in  1971.  He  subsequently  held  post-doctoral  appointments  at 
Battelle  Memorial  Institute  and  Johns  Hopkins  University,  and  most  recently  was 
a visiting  scientist  at  the  Institute  for  Computer  Applications  in  Science  and 
Engineering,  located  at  the  NASA  Langley  Research  Center. 
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WATER  AND  ITS  POSITIVE  ION 


A . Introduction 

The  equilibrium  geometry  of  the  X^A^  state  of  H2O  has  been 
determined  by  infrared  spectroscopy,^  and  various  empirical  deter- 
minations of  the  internal  coordinate  force  constants  have  been 

2 3 4 

derived  from  perturbation  analyses  of  the  spectroscopic  data.  ' ' 

The  deficiencies  in  the  potential  energy  constants  obtained  from 

ah  tKztio  self-consistent-field  (SCF)  calculations  are  apparent 

5 6 

from  the  resultant  computed  vibrational  energies.  ' 

Experimental  equilibrium  geometries  have  been  obtained  for 

22  2 7-10 

the  X B^,  1 Aj^ , and  1 82  states  of  H2O,  but  no  determinations 

of  potential  energy  constants  have  been  made  yet.  Potential  energy 

surfaces  for  the  ground  and  several  excited  states  of  1120''^  have 

been  computed  using  ai  initio  SCF  and  semi-empirical  conf iqurat ion 

interaction  (Cl)  techniques . The  theoretical  rosults^^  and 

12  13 

photoelectron  spectra  ' were  subsequently  combined  to  discuss 

the  predissociation  and  analyze  the  vibrational  structure  of  the 

2 2 
1 B2  state,  and  it  was  concluded  that  the  B2  state  was  predisso- 

9 

dated  twice.  The  position  of  the  appearance  potential  of  the 
+ 2 

OH  ion  being  just  above  the  1 B2  adiabatic  ionization  potential 

is  also  related  to  the  predissociation  of  the  upper  vibrational 

2 2 

levels  of  that  state.  The  coupling  of  the  1 R2  and  1 A^  components 

2 2 
of  the  1 A'  state  adds  to  the  complexity  of  the  B2  surface. 

B.  Neutral  H^O 

One  of  the  authors  (BJR)  has  previously  computed  at  initio 

SCF  and  Cl  potential  energy  surfaces  in  the  neighborhood  of  the 

1 14 

energy  minimum  for  the  X A^  state  of  II2O.  The  energy  surface 
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was  represented  by  a quartic  Taylor  series  expansion,  in  terms  of 
the  internal  coordinates  (AR^,  AR2,  A6),  of  the  form^ 


E = E + 
e 

K^(AR^  + AR2) 

+ K^A0^  + 

2K,AR, AR„ 
6 12 

+ 

2K^{AR^  + AR2)A0  + Kg (AR^  + AR^)  + 

KgA0^ 

+ 

3K^q(AR^  + AR2 

) ARj^AR2  + 

3K^^(AR^  + 

AR2) A0 

+ 

3Ki3(ARi 

+ AR2)A0^ 

+ 

k^^(arJ  + AR2) 

+ K^gAO'^ 

2 

+ AR2)AR^AR2 

+ 

6K^^AR^AR2  + 4K^g(AR^  + 

AR2) A0 

+ 

12K^g(AR^  +AR2 

) AR^AR2A0 

+ 6K2q(Ar2 

+ AR2)A6^ 

+ 12K2j^ARj^AR2A0^  + 4K22<ARj^  + AR2)A9^  (1) 

The  resultant  force  constants  along  with  various  empirically- 
derived  sets  are  given  in  Table  I.  Recently,  in  collaboration  with 
workers  at  other  laboratories,  these  potential  energy  constants  have 
been  used  in  a variation-perturbation  method^  to  compute  vibrational 
energy  levels. The  calculated  and  experimentally-observed^  energy 
levels  are  compared  in  Table  II.  The  importance  of  including  cor- 
relation effects  in  the  wavefunction  when  computing  potential  energy 
surfaces  is  evident  from  the  significant  improvement  in  the  Cl  har- 
monic force  constants  (K^-K^)  and  vibrational  energies  compared  to 
the  SCF  values.  The  Cl  vibrational  energies  show  an  average  error 


of  4%. 
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C.  Positive  Ion 

Ab  initio  SCF  and  multi-configuration  self-consistent  field/ 

configuration  interaction  (MCSCF/CI)  potential  energy  surfaces,  in 

2 

the  neighborhood  of  the  energy  minima,  were  computed  for  the  X 
~ 2 + 

and  B B2  states  of  H2O  . The  energies  were  fit  to  a cubic  Taylor 
series  expansion  (Eq.  (1)  , truncated  after  the  resulting 

equilibrium  geometries  and  force  constants  are  presented  in  Table 
III.  The  computed  equilibrium  geometries  are  seen  to  be  consistent 
with  what  one  would  expect  from  Walsh's  rules.  The  potential  energy 
constants  were  then  input  to  a vibrational  analysis  program^  to  com- 
pute harmonic  frequencies.  Table  IV,  and  second-order  corrected  vi- 
brational energies.  Table  V. 

From  Table  V,  it  is  clear  that  the  computed  vibrational  results 
2 

for  the  B^  state  are  in  good  agreement  with  available  experimental 
7121316 

data.  ' ' ' However,  it  is  seen  that  the  lowest  symmetric  and 

asymmetric  stretching  frequencies  are  stable  to  only  about  10%  as 
one  proceeds  from  the  harmonic  to  cubic  fits,  i.e..  Table  IV  to  Table 
V;  the  lowest  bending  frequency  changes  very  little  which  is  consis- 
tent with  the  fact  that  the  dominant  cubic  term  is  K„ . 

o 

2 

The  results  for  all  the  vibrational  frequencies  of  the  B2 
state  exhibit  significant  changes  as  the  effects  of  cubic  terms  are 
added.  This  is  to  be  expected  at  the  outset,  as  the  anharmonic  nature 
of  the  progression  of  peaks  in  the  photoelectron  spectrum  has  been 
observed. It  is  interesting  to  note  that  the  value  quoted  in  Ref. 


13  for  the  lowest  bending  excitation  agrees  with  that  from  the  cubic 
analysis,  while  the  harmonic  results  agree  with  those  of  Ref.  12. 


r 
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D.  Future  Goals 

Because  of  the  marked  changes  in  the  vibrational  energies  on 

including  cubic  anharmonic  terms  in  the  potential  energy  expansion, 

it  is  necessary  to  investigate  the  effects  of  adding  quartic  and 

higher  terms  in  the  expansion.  In  an  attempt  to  resolve  the  dis- 

2 

puted  assignments  of  the  B2  band  origins  of  the  photoelectron 

spectra  of  H2O,  higher  order  perturbation  theory,  or  a variational 
17  18 

approach  ' will  also  be  considered. 

2 

In  order  to  study  the  predissociation  of  the  B2  state,  it 

will  be  necessary  to  compute  potential  surfaces  for  the  l^A"  and 

2 2 
2 A"  states  as  well  as  additional  regions  of  the  B2  surface.  A 

2 

complete  understanding  of  the  B2  state  will  also  require  an  exam- 
ination of  its  coupling  through  asymmetric  distortion  with  the 
l^A^  state. 
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Table  I.  Equilibrium  Geometry  and  Force  Constants  for  the  X A,  State  of  the  Water  Molecule 
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CARBON  DIOXIDE  AND  ITS  IONS  j 

1 

I 
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A.  Introduction  i 

1 2 

Simple  Walsh-Mulliken  predictions  ' anticipate  much  of  the 
interesting  and  varied  behavior  of  the  CO^ , CO^ / C0~  molecules  in 
a qualitative  sense.  A linear  equilibrium  geometry  is  predicted 

for  ground  state  CO2/  while  many  of  the  low-lying  excited  states  ! 

of  CO2  are  predicted  to  be  bent.  Under  these  conditions,  the  elec- 
tronic spectrum  and  possible  attachment,  detachment,  dissociation, 
etc.,  mechanisms  may  very  widely  with  geometry.  A similar  predic- 
tion is  made  for  CO^.  In  CO^,  however,  there  should  also  bo  many 
low-lying  and  as  yet  unobserved  states  which  arise  from  the  large 
number  of  unoccupied  molecular  orbitals.  The  negative  ion  CO2  is 
predicted  to  have  a bent  equilibrium  geometry.  Again  because  of 
the  large  number  of  unoccupied  valence  molecular  orbitals,  many 
low-lying  excited  states  are  anticipated.  A large  number  of  these 
are  expected  to  be  resonance  states,  inasmuch  as  CO2  will  not  bind 
an  electron  in  many  cases. 

Apart  from  qualitative  considerations  such  as  those  just  out- 
lined, rather  little  information  exists  concerning  the  electronic 

states  of  the  CO2/CO2/CO2  systems.  Definitive  information  consists 

3 4 1 1 

of  Dixon's  assignment  ' of  the  1 B2  1 A^  transition  in  CO2 , 

5 

and  the  four  lowest  vertical  ionization  potentials  of  CO2 • Ryd- 
berg series  in  CO2  have  also  been  accurately  measured,^  but  the 
nature  of  the  Rydberg  orbitals  involved,  and  hence  the  symmetry 
of  the  electron  states,  is  unknown.  The  lower-lying  electronic 
states  of  CO2  are  complicated  by  the  possibility  of  va lence-Rydberg 
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mixing,  and  almost  nothing  is  known  about  these  states.  Even 

the  lowest  valence  type  states  have  not  been  definitely  analy- 

7 8 

zed  experimentally,  and  existing  theoretical  investigations  ' 
of  these  states  have  not  employed  sufficient  accuracy  for  the 
results  to  be  accepted  without  question. 

9 

With  the  exception  of  a semiempirical  study,  little  theo- 
retical work  has  been  done  on  the  electronic  states  of  CO2.  Sub- 
sequent ah  initio  work^^'^^  has  shown  the  semiempirical  results 
to  be  even  qualitatively  incorrect  as  regards  the  stability  of 
ground  state  CO2  relative  to  ground  state  CO2 • The  semiempirical 
model  predicted  that  the  ground  state  potential  energy  surfaces 

of  CO2  and  CO2  do  not  intersect,  whereas  ab  initi'^^'^^  and  exper- 
12 

imental  findings  showed  that  the  surfaces  do  intersect.  As  a 
consequence,  the  assignments  based  on  the  semiempirical  model  are 
open  to  questions. 

B . Linear  and  Bent  Geometry  Correlation  Diagrams 
The  changes  in  the  CO2/CO2/CO2  electronic  states  upon  mole- 
cular bending  were  determined  by  two  sets  of  computations.  First, 
self-consistent-field  (SCF)  calculations  were  carried  out  at  the 
equilibrium  bondlength  (Rj2q=2.1944  bohrs)  of  ground  state  C02(X^"^). 

Second,  the  bondlength  Rqq  was  held  fixed,  and  SCF  calculations  per- 
il 12 

formed  at  a bond  angle  close  to  estimates  ' of  the  equilibrium 

2 

bond  angle  of  ground  state  CO2  (X  A^),  <1  OCO  = 130°.  Results  for 
both  geometries  were  obtained  for  each  of  the  molecules  €02/00^/002. 

Neutral  CO^ 

Figure  1 shows  the  linear  and  bent  geometry  spectra  for  the 

2 • 


lowest-lying  valence  states  of  00 


Most  of  the  excited  states 
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shown  will  have  minimal  energies  for  bent  geometries.  Electronic 
spectra  which  involve  bent  geometry  CO2  molecules  will  be  clearly 
very  different  from  the  vertical  (linear  geometry)  spectrum  of  C02- 
The  bent  molecule  energy  levels  support  Dixon's  assignment  of  B2 
symmetry . ^ A total  of  five  bent  states  (^A^,  ^^2'  ^^2'  ^^2' 
are  bound  relative  to  the  lowest  SCF  asymptote,  + O (^P). 

The  diagram  partially  supports  an  earlier  mechanism^^  for  the  reaction 

0(^P)  + C0(^?:‘^)  ► CO^i^T.'^)  (1) 

in  which  the  surfaces  ^82  and  ^82  are  assumed  to  cross,  and  parti- 
ally supports  a mechanism^  in  which  the  two  surfaces  are  not  assumed 
to  cross.  In  terms  o<  the  former  mechanism,  the  predicted  ^82  surface 
appears  to  overestimate  the  well  depth.  A binding  energy  of  '^-2  cV 
is  predicted  for  ^82  relative  to  0(^P)  + and  this  concurs 

with  both  earlier  kinetic  studies . 

Positive  Ion  CO^ 

Figure  2 shows  the  linear-bent  correlation  diagram  for  CO^ . 

Known  and  unknown  excited  electronic  states  overlap  over  a wide 

energy  range.  The  state  is  substantially  stabilized  upon  ben- 

4 

ding,  where  it  becomes  B^,  and  is  capable  of  predissociating  lower- 

+ 15 

lying  states.  Dissociation  of  CO2  has  been  shown  to  proceed  by 

2 + 

predissociation  of  the  Eg  state 


C02(^E^'^)  + 19.1  eV  > CO^l^Eg'^l 


+ 4 1 + 

o (^S)  + C0(  y ) 


Fluorescence  from  is  not  observed,  and  predissociation  occurs 

2 + 

from  the  ground  vibrational  state  of  Eg  . Our  results  imply  that 

4 2 ^ ‘ 

the  surface  is  capable  of  predissociating  ' in  the*  Franck- 

Condon  region  of  the  ground  vibrational  state.  It  is  furthermore 
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4 + + 

possible  to  correlate  the  Fl^  CO2  and  O , CO  asymptote  symmetries 

with  a C reaction  path.  The  observed  reaction^^ 
s 

C02(^Eg)  + 19.47  eV  CO'^(^Z'^)  + 0(^P)  (3) 

is  suggested  by  our  results  to  proceed  via  a linear  reaction  path, 

4 

again  involving  predissociation  by  the  surface.  The  reaction 

is  reported  to  be  fast.  Our  results  favor  the  asymptotes  CO^ 

2 + 3 

( E ),  0(  P) , since  these  connect  adiabatical ly  to  the  reagent 

+ 4 

states  via  the  CO2 ( surface.  The  set  of  reverse  reactions 

O'^('^S)  + CO(^E''')  --  C0'^(^E''')  + 0(^P)  - .4  eV  (5) 

C0'^(^E'^)  + 0(^P)  ->  00(^2"^)  + O'^('^S)  + .4  eV  (6) 

17-19 

have  been  observed.  According  to  our  results,  the  reagents 

4 " L 

of  Reaction  (!3)  would  approach  on  the  A ( 11^)  surface,  to  which 

they  correlate  in  C symmetry.  Products  would  then  be  formed  by 

4 4 

population  of  the  A' ( H ) surface.  The  exothermic  Reaction (6) 

u 

1 9 

has  been  observed  to  be  fast,  and  could  occur  by  the  reverse  of 
the  process  described  for  Reaction  (5).  However,  there  are  neces- 
sary surface  crossings  in  the  asymptotic  surfaces,  and  for  this 
reason  alone  the  reaction  should  be  fast. 


Negative  Ion  CO^ 

The  linear-bent  correlation  diagram  for  CO2  is  presented  on 

Figure  3.  As  was  found  for  CO^ ^ known  and  unknown  excited  elec- 

4 2 

tronic  states  overlap  over  a wide  energy  range.  The  Hg , 4,  and 

2 2 + 2 + 

Ilg  can  all  be  expected  to  predissociate  the  E and  E states. 


and,  in  some  cases,  possibly  the  Renner  split  component  of 
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2 - 2 

n^.  We  predict  that  C02(X  is  metastable,  in  agreement  with 

other  ah  initio  results , ^ ^ ^ and  with  experiment , but  in  direct 

9 

contradiction  with  semiempirical  findings.  Our  assignments  of 

the  excited  states  also  differ  from  the  semiempirical  assignments. 

2 2 + 

We  find  the  4'g  state  to  lie  above  whereas  the  semiempirical 

study  finds  just  the  opposite.  As  a consequence,  our  correlation 
diagram  supports  different  predissociation  and  reaction  mechanism 
possibilities  from  those  that  could  be  deduced  from  the  semiempiri- 
cal results. 

We  predict  a vertical  electron  affinity  of  ''-4  eV,  which 

20 

agrees  with  experimental  estimates.  The  adiabatic  electron  af- 
finity is  estimated  to  be  '^--l  eV  on  the  basis  of  our  results,  in 
reasonable  agreement  with  experiment . 

The  overall  features  of  the  correlation  diagram  suggest  that 
intersection  of  the  ground  state  CO2/CO2  potential  energy  surfaces 
will  be  similar  to  the  intersection  of  the  isoelectronic  N20/N20~ 
potential  surfaces  which  we  describe  in  detail  in  a later  section 
of  this  report.  We  therefore  predict  that  the  dissociative  attach- 
ment 

C02(X^Eg'*‘)  + e C0(X^E‘'')  + 0"(^P)  (7) 

will  yield  vibrationally  excited  C0(X^E'*^).  This  agrees  with  experi- 

25 

mental  suggestions.  The  associative  detachment 

0“(^P)  + COC^E'*')  •>  C02(^Eg)  + e (8) 

should  also  be  overall  similar  to  the  corresponding  O , N2  detach- 
ment process.  However,  the  electron  affinities  of  CO2  appear  to 

be  significantly  larger  than  those  of  N2O,  and  this  may  account 

27-30 

for  the  fact  that  Reaction  (8)  is  moderately  fast,  whereas 
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the  corresponding  O , N2  detachment  is  slow. 

C . Valence  States,  Mixed  Character  States,  and  Rydberg 
States  in  Neutral  COt 

An  interesting  feature  in  the  electronic  spectra  of  CO2  and 
other  molecules  is  that  the  electronic  states  fall  into  three 
qualitatively  different  classes.  The  first  class  consists  of 
valence  states,  whose  spacial  extents  are  similar  to  that  of  ground 
state  C02(X^Eg)  . The  second  class  consists  of  mixed  character, 

or  large  states,  whose  spacial  extents  are  much  beyond  the  valence 
state  region.  Rydberg  states  make  up  the  third  class,  viz.  Those 
satisfying  the  Rydberg  formula^ 

AE^  = 13.605/[n*(E) (units  are  eV)  (9) 

where 

AE^  = E(Neutral  State)  - E(Parent  Ion  State)  (10) 

and  n*(E)  is  the  effective  quantum  number.  In  the  case  of  mixed 

and  Rydberg  states,  the  character  (spacial  extent,  validity  of 

Equation  (9))  is  determined  by  a single  outer  occupied  orbital. 

We  computed  an  accurate  self-consistent-field  (SCF)  vertical 

spectrum  for  CO2  in  order  to  determine  which  of  the  states  were 

valence,  mixed  character,  or  Rydberg  states.  All  single  electron 

excitations  out  of  the  four  highest  occupied  molecular  orbitals 

(MO's)  into  the  five  unoccupied  MO's  2n  , 5o  , 4o  , 16  and  2 

^ u g u g g 

were  included,  giving  a total  of  60  electronic  states.  The  energy 
levels  are  shown  on  Figure  4.  Orbital  energies,  orbital  second 
moments,  and  the  second  moment  analog  of  Equation  (10), 
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AQj^  = Q(Neutral  State)  - Q(Parent  Ion  State)  (11) 

were  computed.  Effective  quantum  numbers  for  the  mixed  and  Ryd- 
berg states  were  calculated  from  Equation(9),  and  from  the  hydro- 
genic  formula 

n*(r^)  = [0.5bj  + 0.5{b^  + 1.6?^)‘]^  (12) 

with 

b^  = 0.6[lt(^  + l)  - 1/3]  (13) 

The  angular  quantum  numbers,  i,  could  be  determined  without  ambi- 
guity from  the  character  of  the  MO  in  question.  The  SCF  computa- 
tions give  representations  of  five  MO's  (2'"^,  5a^,  4o^,  li*;^  and 
27Tg)  which  are  unoccupied  in  ground  state  C02(X^Zq^).  The  5o^ 
is  always  a large,  or  mixed  character  MO,  and  hence  each  electro- 
nic state  in  which  the  5o^  is  occupied  is  a mixed  character  state. 

7 8 

This  contradicts  earlier,  less  accurate,  ab  initio  findings.  ' 

The  4a^,  16^  and  2tt^  are  always  Rydberg  MO's,  and  thus  electronic 
states  in  which  one  of  these  MO's  is  occupied  are  always  Rydberg 
states.  The  2tt^  MO  is  fundamentally  different,  and  can  have  either 
valence  or  Rydberg  character.  It  has  Rydberg  character  in  the  in- 
states and  ^7!^,  and  otherwise  has  valence  character, 

u g 

The  SCF  and  experimental  Rydberg  series  are  compared  in  Table 
1.  All  experimental  Rydberg  levels  with  n*  < 3 have  been  included, 
since  this  is  the  range  of  n*  we  have  considered  theoretically.  The 
agreement  between  theory  and  experiment  is  very  good,  and  the  sym- 
metries of  the  electronic  states  corresponding  to  the  observed 
series  are  here  definitely  assigned  for  the  first  time.  However, 
none  of  the  SCF  Rydberg  levels  gives  the  experimentally  observed 
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n*  = 2.70.  Preliminary  analysis  suggests  that  the  60^  MO  is  in- 
volved, which  we  did  not  include  in  the  present  work. 

To  a good  approximation,  many  Rydberg  state  characteristics 
depend  only  on  the  Rydberg  MO.  Average  values  for  the  Rydberg  MO's 
are  reported  in  Table  2.  Each  row  of  the  table  contains  the  ave- 
rage value  for  a given  Rydberg  MO.  For  example,  the  second  entry 
in  the  first  row  is  the  average  of  all  the  2it^  MO  orbital  energies 
we  computed.  Deviations  about  the  averages  are  at  worst  a few  per- 
cent, and  usually  smaller.  Note  that  the  two  different  theoretical 

-2 

ways  of  computing  effective  quantum  numbers,  n* (E)  and  n*(r  ),  are 
in  close  agreement. 

D . Vertical  Spectrum  Including  Electron  Corre lat ion  for 

Valence  States  of  the  Configuration  1^^  2’?’^  in  CO-, 

The  configuration  Ift^  can  be  angular  momentum  coupled 

g u 

to  form  the  electronic  states  The  V-state 

is  Rydberg  in  character,  and  will  not  be  considered  here. 

By  allowing  valence  and  Rydberg  character  to  mix,  the  remaining 

states  were  found  to  be  valence  states.  Definitive  experimental 

state  assignments  for  these  valence  states  have  not  been  given, 

7 8 

and  previous  theoretical  treatments  ' of  them  have  been  of  only 
modest  accuracy.  In  the  overall  spirit  of  our  theoretical  project, 
namely,  method  calibration  and  evaluation,  we  have  therefore  per- 
formed several  theoretical  calculations  of  the  electronic  states 
in  question  which  employ  successively  more  accurate  approximations, 
all  of  which  are  more  accurate  than  the  existing  theoretical  data. 
We  first  performed  self-consistent-field  (SCF)  and  multiconfigura- 
tion self-consistent-field  calculations,  followed  by  configuration 
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interaction,  (MCSCF/CI)  with  accurate  (polarized)  expansion  sets. 
These  differed  markedly  from  the  existing  theoretical  spectra, 
implying  that  the  latter  could  only  be  accurate  if  cancellation 
of  expansion  set  and  correlation  energy  errors  occurred.  We  next 
carried  out  equations-of-motion  (EOM)  calculations  with  accurate 
polarized  expansion  sets.  The  EOM  calculations  incorporate  a 
more  complete  treatment  of  the  correlation  energy,  and  provide 
the  most  accurate  theoretical  spectrum  currently  available.  Wo 
find  for  the  EOM  vertical  excitation  energies 


: 

u 

7.41 

eV 

8.13 

eV 

: 

u 

8.45 

eV 

u 

8.51 

eV 

: 

u 

8.51 

eV 

Interestingly  enough,  these  values  are  fortuitously  close  to  the 
existing  limited  accuracy  theoretical  spectra. 

E . Future  Goals 

Potential  energy  surface  priorities  begin  with  the  ground 

31-39 

state  surface  of  neutral  CO2.  More  information  exists  for 

this  potential  surface  than  for  any  others  in  the  molecules  CO2/ 
CO^/  CO2.  We  will  exploit  this  to  simultaneously  assess  our  own 
surfaces  and  several  ways  of  computing  triatomic  force  constants. 
Armed  with  a suitably  calibrated  method,  we  will  then  probe  the 
ground  state  surface  of  neutral  CO2  to  whatever  extent  is  neces- 
sary. Proper  dissociation  will  be  especially  emphasized.  Beyond 
this  we  will  investigate  CO^  potential  surfaces,  which  are  so 
important  in  ionospheric  ion-molecule  reactions.  An  i nvest i qa t i on 


-28- 


of  the  intersection  of  the  CO2/CO2  ground  state  potential  sur- 
faces is  also  envisioned. 

Investigation  of  neutral  CO2  excited  states  will  also  con- 
tinue. Locating  the  lowest  and  states  relative  to  one 

another  appears  to  be  needed  for  understanding  possible  reaction 
mechanisms.  Both  are  optically  forbidden  states,  and  hence  are 
difficult  to  probe  experimentally.  Our  current  theoretical  re- 
sults imply  that  the  states  are  nearly  degenerate.  The  Rydberg 
character  of  the  excited  states  as  a function  of  geometry  is  like- 
wise of  interest.  If  the  Rydberg  character  changes  with  the  geo- 
metry, the  state  in  question  may  not  behave  like  its  parent  ion 
state,  and  hence  its  role  in  chemical  reactions  cannot  be  relia- 
bly predicted  based  on  ion  states. 

Very  little  is  known  about  the  ionization  of  inner  electrons 

4 - 

in  C02<  and  also  of  shake-up  ionization.  The  ion  state  may 

4 

have  predissociation  capabilities  similar  to  the  n ion  state 
that  we  discussed  in  this  report.  The  location  of  the  state, 

and  of  other  unknown  CO^  states,  will  therefore  be  determined. 

There  may  be  low-lying  states  of  CO2  which  lie  in  the  energy 
range  that  we  reported  here.  Using  methods  similar  to  those  we 
describe  for  elsewhere  in  this  document,  a closer  scan  of  pos- 
sible excited  states  will  be  undertaken.  Ultimately,  a re-analysis 
of  the  experiments  as  given  in  Reference  9 is  required  with  reliable 
tth  z t i o curves  . 
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Table  2.  Average  Rydberg  state  expectation  values. 


MO 

AE.^ 

1 

a , b 

-c 

— 2 ^ ^ ^ 
r 

_2C,d 

z 

_2C,d 

X 

AO.® 

1 

n*  (E) 

n*  (r 

2TTg 

1.45 

1.45 

146.3 

62.7 

41.8 

-28.0 

-28.2 

3.07 

3.09 

15 

g 

(d^) 

1.53 

1.53 

120.0 

17.2 

51.4 

46.0 

46.1 

2.98 

2.98 

4o 

u 

2.23 

2.25 

84 . 0 

50.1 

17.0 

-43.8 

-44 . 6 

2.47 

2.51 

2tt 

u 

2.46 

2.49 

66.3 

13.2 

26.6 

17.8 

18.0 

2.35 

2.38 

^Units  are  eV. 

= orbital  energy  of  Rydberg  MO. 

Q 

Rydberg  MO  expectation  value,  origin  at  C and  bonds  along  z-axis. 
d 2 

^Units  are  a^  , a = Bohr  radius, 
o o 

e,,  . . 2.-26 

Units  are  esu-cm  x 10 


■Rydberg  MO  second  moment. 


Fig.  1. 

SCF  State  Correlation  Diagram  for 
Linear  and  Bent  CO2.  The  Energy 
Zero  is  C02(^Sg'''). 


Fig.  2. 

SCF  State  Correlation  Diagram  for 
Linear  and  Bent  COi.  The  Energy 
Zero  is  C02('Lo^)- 
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Fig.  4. 

Vertical  SCF  Spectrum  of  Neutral  CO2  Rela- 
tive to  the  Four  Lowest  States  of  CO^. 
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ozone  AND  ITS  IONS 


A.  Introduction 

Experimentally,  there  are  two  states  connected  by  dipole-allowed 

1 12 

transitions  to  the  ground  A^  state  of  ozone.  ' The  Chappuis 
band  is  assigned  to  the  state,  and  the  Hartley  band  is  assigned 

to  the  ^82  state.  In  addition,  the  Huggins  band  has  been  attributed 
to  the  non-vertical  part  of  the  ^82  ^A^^  transition,  while  the  Wulf 

band  has  been  assigned  to  the  dipole-forbidden  ^A2  ^A^  transition. 

Recent  electron  energy  loss  spectroscopy  has  been  used  to  detect  the 
existence  of  at  least  one  bound  low-lying  triplet  state. ^ The  elec- 
tronic spectra  beyond  the  Hartley  band  has  been  reported,^  including 

4 

the  electron  energy  loss  spectra  for  an  energy  loss  range  of  1-30  eV. 

A variety  of  ab  initio  methods  (self -consistent  field,  generalized 

valence-bond,  and  configuration  interaction)  have  been  used  to  compute 

5 — 8 

the  electronic  spectra  of  ozone.  Three  states  were  found  to  be 

bound  and  a total  of  eight  states  had  excitation  energies  lower  than 
6 eV.® 

The  spectra  of  0^  has  been  studied  by  the  high  resolution  photo- 

9-11 

electron  spectroscopy  of  ozone.  Three  bands,  between  12  and  14 

eV,  with  accompanying  vibrational  structure  were  identified  and  attri- 
buted to  ionization  from  the  6a^,  la2  and  4b2  orbitals.  Two  studies^^ 

2 2 2 

assigned  the  order  of  the  ion  states  as  X A^  , A T\^,  8 82.  whereas  the 
third  study^  assigned  the  ordering  X^A^^,  a'^B2»  B Koopmans'  Theo- 

rem incorrectly  predicts  the  ground  state  of  to  be  ^A2 . with  the 
discrepancy  arising  not  from  differential  orbital  relaxation  effects 
but  from  differences  in  the  correlation  energies  of  the  ionic  states.' 
Multi-configuration  self-consistent  field/configuration  interaction 
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calculations  support  the  assignment  of  the  first  two  excited  states 
2 2 

as  A A2/  B B2  f whereas  the  more  extensive  generalized  valence- 

8 

bond/configuration  interaction  calculations  yields  an  ordering  of 
2 2 

A B2,  B A2.The  photoelectron  spectra  all  show  additional  broad  bands 

at  higher  energies  corresponding  to  ionization  from  inner  orbitals 

followed  by  shake-up  processes,  with  all  three  studies  differing 

on  the  number  and  orbital-state  assignments  of  such  bands.  The  ab 

initio  calculations  also  disagree  with  each  other  (as  well  as  with 

experiment)  on  the  assignment  of  the  higher  states  of  0^. 

Matrix  isolation  spectroscopy  has  been  used  to  study  the  vibra- 

14  15 

tional  and  electronic  spectroscopy  of  the  ozonide  ion.  ' In 

addition,  the  electronic  absorption  and  resonance  Raman  spectra  of 

have  been  studied  by  means  of  y-irradiated  spectra  of  alkali- 
— I — 16,17  ^ J c 


metal  halates. 


Only  one  electronic  transition  is  observed  for 


2 2 14 

and  is  assigned  as  A2  X B^^.  A variety  of  experimental  tech- 

1819  20 

niques  (charge  transfer  reactions,  ' lattice  energy  calculations , 

21  22  23 

collisional  ionization,  drift-tube  photodetachment  ' ) as  well  as 

24 

an  ab  initio  configuration  interaction  calculation  have  been  used 
to  determine  the  electron  affinity  of  ozone,  and  the  consensus  of 
results  indicates  a value  of  2. 0-2.1  eV. 


B.  OZOME 

Ab  initio  self-consistent  field  (SCF)  and  multi-configuration 
self-consistent  field/configuration  interaction  (MCSCF/CI)  calcula- 
tions were  performed  on  the  lowest  singlet  and  triplet  states  of 
at  a single  geometry  corresponding  to  the  experimental  equilibrium 
geometry  of  the  ground  1 state  (Rq_q  = 1.2717A  = 2.403  a^,  0 = 
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2 5 2 6 

116.78°)  . Based  on  earlier  SCF  calculations,  the  manifold  of 

states  of  a given  symmetry  which  result  from  various  occupancies 

of  the  6a^,  4b2/  and  la2  valence  orbitals  are  all  quite  close 

in  energy.  Thus  a priori  selection  of  the  dominant  Hartree-Fock 

configuration  for  a given  state  is  not  possible,  and  it  is  necessary 

to  compute  energies  for  several  configurations  of  the  same  symmetry 

in  order  to  determine  the  lowest  one. 

The  computed  SCF  and  MCSCF/CI  vertical  excitation  energies 

for  the  valence  states  of  are  given  in  Table  I.  An  examination 

3 

of  the  SCF  results  reveals  that  the  B2  state  is  incorrectly  predicted 
to  be  the  ground  state  at  this  level  of  computation;  this  incorrect 

g 

state  ordering  was  also  found  in  other  theoretical  calculations. 

Moreover,  there  is  a reordering  of  some  of  the  other  excited  states. 

The  including  of  correlation  energy  results  in  a correct  ordering  of 
1 3 

the  and  B2  states.  The  present  MCSCF/CI  excitation  energies 

are  within  the  reported  experimental  range  of  energies  for  the  ob- 
served transitions.  It  should  be  noted  that  the  order  and  position 
of  the  electronic  states  of  0^  can  be  expected  to  vary  siani f icant ly 
with  geometry,  even  to  the  extent  of  inverting  the  ground  and  excitc<) 

4-4-  8 

states . 


C.  POSITIVE  ION 

SCF  and  MCSCF/CI  calculations  were  performed  on  the  lowest  dou- 
blet and  quartet  states  of  O^  at  the  experimental  equilibrium  geome- 
try of  the  X^A^  state  of  O^-  The  results  of  the  vertical  excitation 
energies  are  given  in  Table  II.  On  theoretical  grounds  we  expect 
0^  to  have  an  energy  surface  behavior  comparable  to  NO2.  Although 
Koopmans'  Theorem  predicts  the  ground  state  of  0“^  to  be  '^^2' 
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experimental  ground  state  of  the  isoelectronic  NO2  molecule  is 
2 1 

A^.  Once  again  the  SCF  calculations  do  not  yield  the  correct 

12  13 

ground,  as  also  found  in  other  theoretical  works.  ' Including 

correlation  energy  results  in  the  correct  ground  state  and  shov/s 

2 

the  first  excited  state  to  be  62- 

The  computed  three  lowest  ionization  potentials  of  are 
given  in  Table  III,  and  are  compared  with  experimental  results. 

The  MCSCF/CI  values  are  all  in  error  by  1.5-1. 8 eV.  We  presently 
attribute  the  large  errors  in  our  calculated  ionization  potentials 
to  restrictions  in  the  present  level  of  our  com.putat iona  1 method. 

D.  NEGATIVE  ION  0~ 

SCF  and  MCSCF/CI  calculations  were  performed  on  the  lowest  dou- 
blet and  quartet  states  of  0^  at  the  experimental  eauilibrium  geo- 

1 14 

metry  of  the  X A^  state  of  O^.  Previous  experimental  and  theore- 

24  2 

tical  studies  estimate  the  equilibrium  geometry  of  the  X state 

of  to  be  near  that  of  the  ground  state  of  the  neutral  molecule. 
The  computed  vertical  excitation  energies  are  given  in  Table  IV.  In 
contrast  to  what  was  found  for  the  and  systems,  SCF  calcula- 
tions are  sufficient  for  a correct  prediction  of  the  ground  state 
of  O^.  The  results  for  the  quartet  vertical  excitation  enerqies 
must  be  viewed  with  caution  since  these  states  are  unstable  in  the 

SCF  approximation.  The  MCSCF/CI  calculation  is  in  error  by  1 eV 

2 2 

with  the  one  experimentally  observed  transition,  ^2  * ^ 

present  results  would  seem  to  better  support  an  assignment  of 
2 2 

Aj^  • X Bj^,  but  deficiencies  in  the  computations  preclude  any 
reassignments  at  the  present  time.  It  should  also  be  noted  that  the 
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r 

I 

electronic  spectra  of  0^  were  studied  in  an  arqon  matrix,  not  in 
' the  gas  phase;  thus  the  value  of  the  measured  transition  energy 

may  be  significantly  affected  by  environmental  interactions,  '^he 

[. 

1 computed  SCF  electron  affinity  exceeds  the  experimental  values  by 

' 0 . 3 eV,  while  the  use  of  correlated  wavef unctions , in  their  pre- 

i 

j sent  level,  fails  to  bind  the  electron  to  ozone. 

i 

E.  FUTURE  GOALS 

Future  work  will  deal  with  the  extension  of  present  levels  and 
i types  of  computational  procedures  in  order  to  more  accurately  cal- 

I culate  energy  quantities,  especially  energy  differences  between 

systems  having  different  numbers  of  electrons  (e.g.,  ionization 
potentials,  electron  affinity) . A primary  objective  of  this  work 
will  be  the  assignment  of  the  high  energy  bands  in  the  photoelectron 
spectra  of  O^,  as  well  as  the  prediction  of  some  yet  unobserved  tran- 
I sitions  in  the  0^  ozonide  ion. 

I 

I We  also  propose  to  compute  portions  of  the  potential  energy  sur- 

face of  the  X^Aj^  state  of  ozone,  with  the  principal  emphasis  in  the 

well  region  to  compute  bound  vibrational  levels  and  spectroscopic 

3 - 3 27-21 

constants.  The  vibrational  deactivation  of  ( y^)  by  0(  P) , 
and  the  three  body  recombination  process, 

0(^P)  + 02(^i;g)  + M ->■  03(X^Aj^)  + 

«■  are  additional  topics  of  interest  involving  the  ground  state  poten- 


tial surface. 
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Table  II. 

Vertical  excitation 
0 = 116. 8”.^ 

energies  for 

Vo  - "•‘">5  “o- 

State 

SCF 

MCSCF/CI 

3.8 

0.0 

0.7 

0 . 5 

5.6 

3.2 

3.8 

0.1 

l\ 

1.7 

1.9 

0.0 

2.7 

l\ 

4.9 

4.9 

1^2 

1.0 

1.4 

a 


All  energies  in  electron  volts. 
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Table  IV.  Vertical  excitation  energies  for  0^,  Rq_q  = 2.403  a^, 

e = 116.8“.® 


State 

SCF 

MCSCF/CI 

Expt . 

0.0(2. 4) 

0.0(-0.5) 

0.0(1.47- 

l^Ai 

2.6 

2.7 

1^A2 

4.2 

3.7 

2.8 

l'«2 

2.9 

2.9 

l^Ai 

10.1 

10.8 

i'*A2 

9.0 

8.3 

l^Bi 

9.3 

8.3 

1^2 

7.3 

8.3 

All  energies  in  electron  volts. 

^M.  E.  Jacox  and  D.  E.  Milligan,  J.  Mol.  Spectry.  43,  148  (1972). 
'Number  in  parentheses  is  the  electron  affinity. 

^J.  Berkowitz,  W.  A.  Chupka,  and  D.  Gutman,  J.  Chem.  Phys.  55,  2733 
(1971). 

'E.  E.  Ferguson,  Can.  J.  Chem.  47,  1815  (1969). 

^R.  H.  Wood  and  L.  A.  D'Orazio,  J.  Phys.  Chem.  69,  2562  (1965). 


^E.  W.  Rothe,  S.  Y.  Tang,  and  G.  P.  Reck,  J.  Chem.  Phys.  62,  3829 
(1975). 

*^R.  Byerly,  Jr.,  and  E.  C.  Beaty,  J.  Geophys.  Res.  76,  4596  (1971). 

^S.  F.  Wong,  T.  V.  Vorburger,  and  S.  B.  Woo,  Phys.  Rev.  A5,  2598  (1972). 

^M.M.  Heaton,  A.  Pipano, and  J.J.  Kaufman,  Int.  J.  Quantum  Chem.  56,181(1972! 
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NITROUS  OXIDE  AND  ITS  NEGATIVE  ION 
A.  Introduction 

The  N2O,  N2O  molecules  are  isoelectronic  with  the  systems  | 

12 

CO21  CO2#  respectively,  and  the  Walsh-Mulliken  rules  ' predict  1 

that  linear  and  bent  geometries  play  important  roles  in  N2O,  N2O  | 

chemistry.  The  ground  state  N2O/N2O  prediction  is  that  the  mini- 
mum geometry  of  N2O  is  linear,  and  the  minimum  geometry  of  ground 
state  N^O  is  bent.  As  was  the  case  in  the  CO2 , CO2  molecules, 
the  potential  energy  surfaces  intersect,  with  the  ground  state 
potential  energy  surface  of  N2O  lying  lowest  for  geometries  near 

the  linear  equilibrium  geometry  of  N2O,  and  with  the  ground  state  : 

potential  energy  surface  of  N2O  lying  lowest  near  the  bent  equi- 
librium geometry  of  N2O  . The  details  of  these  surfaces  and  their 
intersection  bear  directly  on  the  possible  ionospheric  associative 
detachment  reaction 

0”(^P)  + N2(X^Eg)  N20(X^E'*')  + e + 0.21  eV  (1) 

which  has  not  been  observed  in  laboratory  experiments  conducted 
at  300K.  A low  thermal  rate  constant  has  been  assigned  to  Reaction 
(1)  by  several  experimenters,  and  conflicting  reports  exist'^' 
concerning  reactivity  above  0.3  eV  relative  collision  energy  (ground 
vibrational  state  N2) . 

The  low  thermal  rate  constant  for  the  exothermic  Reaction  (1) 
suggests  that  substantial  barriers  are  involved  in  the  intersection 
of  the  N2O  and  N2O  potential  energy  surfaces,  and  that  vibrational ly 
excited  reagent  N2  may  appreciably  enhance  the  reaction.  Since  vib- 
rationally  excited  N2  is  common  in  the  ionosphere,  the  characteriza- 
tion of  the  N2O  and  N2O  potential  surfaces,  needed  to  predict  the 

_J 


-49- 


1 

) 

effects  of  vibrationally  hot  reagent  N2  on  the  rate  of  Reaction  j 

{!),  warrants  detailed  investigation . | 

B . The  Ground  State  N^O  and  N^O  Potential  Energy  Surfaces 
Thermodynamic  cycles  can  be  constructed  which  express  the 
adiabatic  electron  affinity  of  N2O,  EA(N20),  and  the  dissociation 
energy  of  the  N-NO  bond  in  N2O  , D(N-NO  ) , in  terms  of  known  ther-  1 

modynamic  quantities:  ] 


EA(N20) 

= EA(0)  + D(N2-0  ) - D(N2-0) 

(2a) 

D (N-NO“ 

) = EA(N20)  + D(N-NO)  - FA (NO) 

(2b) 

The  values  used  are  shown  in  Table  1,  together  with  the  resulting 
EA(N20)  and  D(N-NO  ) values. 

Ab  initio  self-consistent-field  (SCF)  and  multiconfiguration 
self-consistent-field,  followed  by  configuration  interaction  (MCSCF/ 
Cl) , wavefunction  calculations  were  performed  to  determine  the  pot- 
ential energy  surfaces  for  N2O  and  N2O  . Equilibrium  geometries 
and  harmonic  force  constants  and  vibrational  frequencies  are  shown 
in  Table  2.  It  follows  from  these  results  that  the  shapes  of  the 
N2O  and  N2O  potential  energy  surfaces  should  be  reliably  determined 
near  their  respective  minima. 

Ab  initio  bending  and  stretching  potentials  can  be  used  to  esti- 
mate the  energy  required  to  change  the  geometry  of  N20~  from  its 
equilibrium  value  to  that  of  equilibrium  N2O.  The  computed  values 
are  listed  in  Table  3.  The  predicted  barrier,  2.5  eV,  is  signifi- 
cantly higher  than  the  exothermicity  of  Reaction  (1),  in  support  of 


a low  thermal  rate  constant. 
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The  intersection  locus  of  the  ground  state  N2O  and  N2O  poten- 
tial energy  surfaces  can  be  determined  by  fitting  the  ab  initio 
N2O  and  N2O  potential  energy  points  with  a functional  form.  It 
was  found  that  a two-dimensional  anharmonic  Morse  potential  form 
was  sufficient,  in  which  one  bondlength  is  treated  parametrically. 
The  anharmonicity  is  essential  for  close  reproduction  of  the  ab 
initio  N^o  vertical  electron  affinity.  The  minimum  energy  of 
intersection  between  the  N2O  and  N2O  potential  surfaces  is  found 
to  be  ^^^0.66  eV  above  the  potential  energy  of  equilibrium  geometry 
N2O.  The  geometry  at  the  minimum  energy  of  intersection  is 


O 


(3a) 

■^N  ■ 

(3b) 

<)(NNO)  = 155° 

(3c) 

The  associative  detacliment  threshold  for  Reaction  (l)is  obtained 
from  the  minimum  energy  intersection  by  adding  the  zero  point  energy 
of  the  minimum  intersection  energy  complex  and  subtracting  1)  the 
asymptotic  N2  + O energy  obtained  from  the  anharmonic  Morse  form, 
and  2)  the  zero  point  energy  of  N2*  The  predicted  threshold  for 
Reaction  (1)  is 

E ( Rxn  1 ) = 0 . 2 1 eV  { 4 ) 

d 

This  activation  energy  is  less  than  the  energy  necessary  to  excite 
reagent  N2  to  the  first  excited  vibrational  state  (0.29  eV) . Fur- 

O 

thermore,  the  of  Equation  (3b)  is  greater  by  0.04  A than  the 

O 

classical  anharmonic  turning  point  for  vibrat ionally  cold  N2(1.146A), 
and  almost  equal  to  the  outer  turning  point  for  N2  in  its  first 

O 

excited  vibrational  state  (1.185A).  Thus,  the  rate  of  Reaction  (1) 
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will  probably  not  be  enhanced  by  translational  excitation  of  vib- 
rationally  cold  N2/  but  may  be  strongly  enhanced  by  vibrational 
excitation  of  reagent  N2»  with  or  without  translational  excitation. 

The  N2O  and  N2O  potential  energy  surface  intersection  also 
provides  information  for  the  dissociative  attachment  threshold 

N20{X^E'*')  + e ->■  N2(X^Eg)  + o“(^P)  - .21  eV  (5) 

The  electron  detachment  threshold  for  Reaction  (5)  is  obtained 
from  the  minimum  energy  intersection  by  adding  the  zero  point  energy 
of  the  minimum  intersection  energy  complex,  and  subtracting  1)  the 
N2O  equilibrium  potential  energy,  and  2)  the  zero  point  energy  of 
N2O.  The  result  is 

Ea(Rxn(5) ) = 0.40  eV  (6) 

0 

which  is  in  good  agreement  with  the  experimental  estimate  of  0.45  eV. 
Since  the  minimum  intersection  occurs  with  extended  bondlengths  rela- 
tive to  equilibrium  N2O,  and  at  an  angle  about  25°  away  from  lin- 
earity, it  is  very  likely  that  the  electron  detachment  described  by 
Reaction  (4)  will  be  facilitated  by  excitation  of  the  N2O  bending 
and  symmetric  stretch  modes. 

The  vertical  electron  affinity,  or  resonance  energy,  for  N20(^ 
can  be  reliably  obtained  from  the  previously  computed  adiabatic 
electron  affinity  by  adding  the  zero  point  energy  of  N2O  , and  sub- 

- 2 

tracting  1)  the  zero  point  energy  of  N2O  ( II)  , which  is  estimated 
to  be  the  same  as  the  N2  zero  point  energy,  2)  the  barrier  to  ben- 
ding, and  3)  the  barrier  to  contraction  of  the  bondlengths.  A 

vertical  electron  affinity  of  -2.23  eV  results,  which  is  in  good 

q 

agreement  with  broad  experimental  electron  scattering  peaks. 


I 
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C . Conclusions 

The  ground  state  potential  energy  surface  of  N2O  is  stable 
in  its  equilibrium  region  to  both  associative  detachment  (Reaction 
(1))  and  dissociation.  The  features  of  the  ground  state  N2O/N2O 
potential  surfaces  described  in  this  report  show  why  this  is  so. 
Associative  detachment  is  immeasurably  slow  for  reagent  thermal 
N2  because  there  is  a barrier  of  vO.2  eV  which  the  reagents  must 
overcome  to  reach  the  minimiom  intersection  energy  of  the  N2O/N2O 
surfaces.  Moreover,  the  minimum  intersection  occurs  at  bond- 

lengths  greater  than  the  equilibrium  N2  bondlength,  which  suggests 
that  vibrationally  hot  N2  will  be  needed  to  effectively  enhance 
reaction . 

The  dissociative  attachment  process  (Reaction  (5))  is  likewise 
described  by  the  M2O/N2O  potential  surfaces.  The  N2O  molecule 
must  be  vibrationally  excited  in  order  to  reach  the  minimum  energy 
of  intersection  of  the  N2O/N2O  surfaces.  This  requires  at  least 
the  absorption  of  several  bending  vibrational  quanta.  Here,  how- 
ever, the  threshold  energy  (0.4  eV)  is  sufficient  to  populate  these 
levels.  Since  the  minimum  intersection  energy  involves  slightly 
longer  R^^^  and  R^q  than  occur  in  N2O,  it  is  likely  that  the  exci- 
tation of  stretching  vibrational  quanta  will  further  enhance  reac- 
tion . 

Our  results  also  provide  a vertical  electron  affinity  of  -2.23 
eV  for  N20(^1^).  This  value  is  thought  to  be  reliable,  since  it 
was  obtained  with  a method  which  does  not  involve  a nonvariat ional 


determination  directly. 
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D.  Future  Goals 

There  is  evidence  from  experimental and  theoretical ^ ^ 
studies  that  a major  source  of  N2  vibrational  temperature  in  the 
upper  atmosphere  is  its  E -»•  V,R,T  energy  transfer  reaction  with 
0(^D) . The  vibrational  temperature  of  N2  is  important  because 
it  strongly  influences  important  atmospheric  characteristics  such 
as  electron  density  and  temperature,  and  the  quantity  of  infrared 
radiation  from  CO2,  H2O,  etc.  The  intersection  of  excited  N2O 
triplet  surfaces  with  the  ground  state  singlet  N2O  surface  will 
be  determined,  since  these  are  central  to  understanding  the  vib- 
rational deactivation  of  N2  by  0(^D) . Preliminary  SCF  computa- 
tions in  our  laboratory  indicate  that  correlated  MCSCF/CI  wave- 
functions  are  needed  for  these  surfaces. 

The  reaction  to  form  NO^{X^Z^)  + N(^S)  from  O^(^S)  + N2(X^X^) 

14 

determines  the  free  electron  depletion  rate  in  the  ionosphere. 

The  understanding  of  this  reaction  involves  the  intersection  of 
excited  N20^  quartet  surfaces  with  the  ground  state  N20^  doublet 
surface.  We  shall  investigate  these  in  detail,  since  there  may 
be  reaction  barriers  which  affect  the  accessible  paths. 
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Table  1.  Bond  dissociation  energies  and  electron  affinities  for  the  N2O/N2O  systems. 
Units  are  eV. 
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Table  2.  Ab  initio  potential  energy  surface  characteristics  for 
N2O  and  N2O  . 


Parameter 

N.,0 

N^O 

SCF 

MCSCF/CI 

MCSCF/CI 

Expt . 

1.204 

1.222 

1 .147 

1.128^ 

1 . 376 

1 . 375 

1.240 

1.184^ 

<)NNO(deg) 

123.6 

132.68 

180 

180^ 

kNN(md/A) 

11.54 

11.49 

17 .21 

17.88^ 

kNo(md/A) 

3.93 

3.83 

9.99 

11.39^ 

0.645 

0.643 

0.41 

0.49^ 

/ ~1\  c 

VNo(cm  ) 

945 

912 

1223 

1285^ 

535 

555 

524 

589^ 

1660 

1666 

2183 

2224^ 

0.196 

0.195 

0.245 

0.254^ 

^G.  Herzberg,  "Electronic  Spectra  and  Electronic  Structure  of 
Polyatomic  Molecules"  (Van  Nostrand  Reinhold,  New  York,  1966). 

^E . B.  Wilson,  J.  C.  Decius,  and  p,  C.  Goss,  "Molecular  Vibrations 
(McGraw-Hill,  New  York,  1955). 

^Atomic  masses  are  mj^  = 14.00751,  m^  = 16.0000 

d _ . . 

Zero  point  energy. 

\ 


I 

1 

! 
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Table  3.  Bending  and  stretching  barriers  for  N2O 


SCF 

MCSCF/CI 

AEj^(eV)^ 

1.41 

1.10 

AEs(eV)^ 

-1.00 

0 

AE(eV) ^ 

2.41 

2.50 

^Energy  rise  upon  bending  (NNO)  from  the  bent  equilibrium  value 
to  180°.  Bondlengths  are  equal  to  their  equilibrium  values. 

^Energy  change  upon  stretching  the  NN  and  NO  bondlengths  from  their 
equilibrium  N2O  values  to  their  equilibrium  N2O  values. 

^AE  = AE,  - AE^ 

D S 
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lONS  OF  NITROGEN  DIOXIDE 

A.  Introduction 

High  resolution  photoelectron  spectra  for  NO2  are  known  for 

12  + 

energies  up  to  27.5  eV.  ' The  NO2  ion  peaks  have  been  assigned 

on  the  basis  of  ground  state  NO^  calculations  and  by  analogy  with 

2-4 

the  isoelectronic  CO2  molecule.  No  definitive  experimental 

assignment  of  the  first  ionization  potential  has  been  made,  owing 
to  the  great  difference  between  the  neutral  (bent)  and  ion  (linear) 
equilibrium  geometries. 

Semiempirical  calculations  on  NO2  are  in  essential  agree- 

ment with  a number  of  experiments  dealing  with  absorption  spectra 
in  NaM02  or  phosphorescence  in  aqueous  solutions.  Polarization 

studies  have  assigned  the  electric  dipole  allowed  state  in  the 

7 8 3 

region  just  under  3 eV.  ' The  state  has  been  assigned  to 

9-11  1 

v2.4  eV.  The  state  was  calculated  to  be  v4 . 5 eV  semiempi- 

rically,^  and  on  this  basis  spectral  assignments  made . The 
^82  state  is  dipole  allowed,  and  lies  about  6 eV  above 

B.  Neutral  NO^ 

The  most  accurate  theoretical  results  available  for  NO2  were 

14 

previously  calculated  in  this  laboratory.  We  use  these  to  ev^al- 
uate  our  present  efforts  for  the  singly  charged  ions.  The  vertical 
NO2  spectrum  we  compute  agrees  with  the  accurate  results  to  within 
'vO.B  eV.  To  the  extent  that  NO2  is  representative  of  NO^  and  NO^  » 
a comparable  accuracy  of  the  ionic  spectra  is  expected. 
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C.  Positive  Ion  NO^ 

Vertical  excitation  energies  for  NO2  = 2.25  bohrs,  <) 

ONO  = 134°)  were  computed  with  single-  and  multiconfiguration 
self-consistent-field  methods  (SCF  and  MCSCF  methods) . Config- 
uration interaction  was  then  carried  out  using  the  MCSCF  results 
(MCSCF/CI) . The  computed  and  experimental  spectra  are  compared 
in  Table  1.  The  SCF  results  predict  the  wrong  first  ionization 
potential  (^^2  instead  of  ^A^) , but  otherwise  are  in  reasonable 
agreement  with  the  MCSCF/CI  spectrum.  The  MCSCF/CI  spectrum  is 
somewhat  different  from  the  experimental  spectrum,  and  this  may 
in  part  be  due  to  the  difference  between  theoretically  computed 
vertical  ionizations  and  those  obtained  by  subtracting  experi- 
mental photoelectron  spectral  peaks.  This  is  probable  for  mole- 
cules like  N02/N02f  where  ground  state  neutral  and  ion  equilibrium 
geometries  are  very  different.  Because  of  these  geometry  diffe- 
rences, vibrationally  excited  levels  can  be  expected  to  be  obser- 
ved in  the  experiments,  and  these  are  not  included  in  theoretical 
calculations  such  as  ours. 

A linear  to  bent  geometry  NO^  correlation  diagram  is  drawn 
on  Figure  1.  The  energies  represented  by  solid  lines  are  SCF  ener- 
gies. The  dashed  lines  represent  MCSCF  variational  estimates  of 
states  which  collapse  variationally  in  the  SCF  model.  The  dotted 
line  shown  for  (^82)  should  not  be  taken  too  seriously,  since  it 
is  a guess  made  on  the  basis  of  the  Walsh-Mulliken  rules, 
and  was  not  observed  in  our  computed  results.  The  linear  NO^  SCF 
spectrum  appears  to  be  correctly  ordered,  and  clearly  shows  the 
greater  stability  of  the  linear  ground  state.  However,  the  beha- 
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vior  of  the  states  introduced  by  bending  is  not  like  that  obser- 
ved in  C02>  this  is  possibly  significant  because  in  the  past 

orderings,  etc.,  have  been  rationalized  by  assuming  CO2  is  a rea- 
sonable model  for  NO^- 

D.  Negative  Ion  NO^ 

Computed  SCF  vertical  excitation  energies  (Rj,jq  = 2.25  bohrs, 
ONO  = 134°)  and  experimental  excitation  energies  for  NO~  are 
shown  in  Table  2.  Qualitative  agreement  is  observed  between  the 
SCF  and  experimental  spectra,  but  quantitative  differences  are 
obvious.  Reasons  for  the  quantitative  differences  may  in  part 
be  due  to  differences  between  the  vertical  geometry  and  the  experi- 
mental geometry.  For  example,  the  <)  ONO  for  the  minimum 

17  18 

geometry  may  be  'v,  116°  ' as  opposed  to  the  vertical  value 
<)  ONO  = 134°.  Another  possible  reason  is  that  the  SCF  method  is 
not  totally  appropriate  for  the  kinds  of  vertical  electron  attach- 
ments that  occur  in  forming  NO2  from  NO2 . 

An  adiabatic  correlation  diagram  for  NO^  is  drawn  on  Figure  2. 
The  C2^  energy  levels  are  those  reported  in  Table  2.  The  diagram 
on  Figure  2 is  the  first  for  NO2  that  is  based  on  ah  initio  enerqy 
levels . 

According  to  the  correlation  diagram,  electron  attachment  in 

3 1 

the  range  %0-l  eV  is  associated  with  only  two  states,  and  X A^. 

Structure  which  is  observed  below  1 eV  should  be  associated  with 

19-21 

the  ground  state,  in  agreement  with  experimental  conclusions. 

Dissociative  attachment  of  electrons  has  produced  each  of  the  spe- 
- - - 22-24 

cies  O , O2  and  NO  . The  appearance  potentials  for  each  are 
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about  the  same  as  their  thermodynamic  thresholds.  The  adiabatic 

paths  on  Figure  2 likewise  suggest  no  apparent  barriers.  The  pro- 

22 

duction  of  O is  observed  to  be  largest,  and  this  agrees  with 
the  fact  that  several  adiabatic  paths  on  Figure  2 provide  routes 
to  the  lowest  energy  asymptote,  which  produces  O . However,  non- 
adiabatic  interactions  must  occur  for  the  case  of  NO  production. 

E . Futxire  Goals 

The  SCF  linear  to  bent  geometry  correlation  diagrams  for  NO^ 
and  CO2  which  we  have  computed  in  the  present  work  suggest  that 
the  isoelectronic  CO2  and  NO^  molecules  behave  differently  in  a 
qualitative  sense.  We  shall  determine  whether  this  difference  pre- 
vails in  more  accurate  descriptions.  Several  of  the  states  which 
enter  in  NO^  Walsh  rules  will  also  be  computed  in  order  to  furnish 
a completely  ab  initio  correlation  diagram.  The  extended  and  more 
accurate  descriptions  of  NO^  will  also  be  needed  to  further  under- 
stand the  difference  between  theory  and  experiment,  if  indeed  this 
difference  persists. 

Negative  ion  NO2  energy  levels  which  include  electron  corre- 
lation, as  well  as  different  bending  angles,  are  likewise  needed  to 
clarify  the  differences  between  theory  and  experiment.  In  this 
respect  NO2  is  a good  candidate  system,  since  it  is  variationally 
stable  relative  to  NO2  in  many  states,  and  since,  experimentally, 
optically  allowed  transitions  are  low-lying  in  the  NO2  spectrum. 
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Table  1.  Vertical  excitation  energies  for  NO^/  = 2.24  bohrs, 

<)  ONO  = 134°.  Units  are  eV. 


Ion  State 

SCF 

MCSCF/CI 

Expt . ^ 

0.90 

0.00 

0.00 

'®2 

0.00 

2 . 56 

1.78 

3a, 

1.34 

2.65 

2.37 

'^2 

1.78 

3.55 

2.83 

4.13 

3.84 

3.28 

7.28 

6.83 

6.22 

8.27 

7.70 

6.41 

9.28 

8.78 

- 

^C.  R.  Brundle,  D. 
and  D.  G.  Streets 

Neumann , 

, J.  Chem. 

W.  C.  Price,  D. 
Phys . 53,  705 

Evans,  A.  W.  Potts, 
(1970)  . 

-65- 


Table  2. 

Excitation  energies 
134°.  Units  are  eV. 

for  NO2, 

- 2.2  5 bohrs,  <)  ONO  = 

Ion  State 

SCF^ 

Expt 

'^1 

0.00 

0.98 

2.47 

'^2 

4.12 

'^2 

5.60 

"^2 

5.92 

. 5^ 

"^2 

8.26 

9.62 

^Vertical 

^Adiabatic 

^H.  J.  Maria,  A.  Wahlborg,  and  S.  P.  McGlynn,  J.  Chem.  Phys.  49, 
4925  (1968). 

R.  M.  Hochstrasser  and  A.  P.  Marchetti,  J.  Chem.  Phys.  50,  1727 
(1969). 
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c. 

Allen  and 
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Dixon, 

Trans.  Faraday 

Soc.  65,  1168  (1969) 

J. 

w. 

Sidman,  J. 

Amer . 

Chem. 

Soc.  79,  2669 

(1957)  . 

J. 

w. 

Sidman,  J. 

Amer . 

Chem . 

Soc.  79,  2675 

(1957)  . 
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McEwen,  J. 
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Phys . 
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^L.  E.  Harris,  H.  J.  Maria,  and  S.  P.  McGlynn,  Czech.  J.  Phys.  B20, 
1007  (1970) . 


-66- 
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Fig.  I. 

State  Correlation  Diagram  for  Linear  and 
Bent  NO'*’.  The  Solid  Lines  are  SCF  Energy 
Levels.  The  Dashed  Lines  are  Variational 
MCSCF  Values  for  Higher  Roots  of  a Given 
Symmetry.  The  Dotted  Value  is  an  Estimate 
Based  on  the  Wal sh-Mul I i ken  Rules,  and  was 
not  Computed  ab  initio.  The  Energy  Zero 
is  NO^C^Lg''')  as  Given  by  the  SCF  Model. 


Fig.  2. 

State  Correlation  Diagram  for  NO2.  The 
Molecular  Ion  Energies  are  SCF  Values.  The 
Asymptote  Energies  were  taken  from  "Poten- 
tial Energy  Surfaces  for  Air  Triatomics," 
ii.  Krauss,  D.  G.  Hopper,  P.  J.  Fortune, 

A.  C.  Wahl,  and  T.  0.  Tiernan,  ARL  TR  75- 
0202,  Vol . I.  The  Energy  Zero  is  N02(^A]) 
as  Given  by  the  SCF  Model. 
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